Anisotropic J/ty suppression in nuclear collisions 
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The nuclear overlap zone in non-central relativistic heavy ion collisions is azimuthally very 
asymmetric. By varying the angle between the axes of deformation and the transverse direction of 
the pair momenta, the suppression of J/$ and \£ f/ will oscillate in a characteristic way. Whereas 
the average suppression is mostly sensitive to the early and high density stages of the collision, the 
amplitude is more sensitive to the late stages. This effect provides additional information on the J/^ 
suppression mechanisms such as direct absorption on participating nucleons, comover absorption or 
formation of a quark-gluon plasma. The behavior of the average J/^ suppression and its amplitude 
with centrality of the collisions is discussed for SPS, RHIC and LHC energies with and without a 
phase transition. 

I. INTRODUCTION discussed in Sec. V, and finally a summary is given. 



The observed suppression of J/^f and in nuclear 
collisions at SPS energies [0 has created much inter- 
est, in particular whether the "anomalous" suppression 
in central Pb + Pb collisions can be attributed to the 
formation of a quark-gluon plasma (2|^]. It has been 
found (see, e.g., Q) that in p+A reactions the suppres- 
sion can be understood as Glauber absorption of Char- 
monium states ip = J/^,^f',X; on participant nucleons 
with a\s/N — 7.3 mb. At SPS energies also the J/flj sup- 
pression in S + A and peripheral Pb+Pb collisions can be 
explained reasonably well but the model fails for central 
Pb+Pb reactions. Additional break up on comovers can, 
however, also account for the observed strong suppres- 
sion H as is supported by microscopic approaches Hf7|. 
Understanding the early state dynamics in ultrarelativis- 
tic heavy ion reactions is one of the major tasks in view 
of the upcoming experiments at RHIC and LHC. 

We will here suggest a new method to explore the J/*S> 
suppression mechanism which exploits the anisotropy 
of the overlap zone in semi-central collisions. We will 
demonstrate how the azimuthal distributions of charmo- 
nia depends on the initial asymmetry and can give ad- 
ditional information on formation/dissociation processes, 
in particular when and where they occur. We find that 
the azimuthal asymmetry is ^5-10% for semi-central 
Pb+Pb collisions at SPS energies and more than doubles 
at RHIC and LHC energies. This is to be compared to 
the present experimental uncertainty of 1-2% in the J/ty 
minimum bias data of NA50 [Q. It is, however, less than 
the typical elliptic flow 2v 2 ~ 20% at SPS @. Measuring 
the reaction plane in coincidence with J/'f's is therefore 
a promising method which provides new insight in J/\f' 
absorption mechanisms and possible phase transitions. 

The two important mechanisms for J/\l/ suppression, 
comover and Glauber absortion, are discussed in Sec. 
II and III respectively with emphasis on the azimuthal 
asymmetries. Result for SPS, RHIC and LHC energies 
are presented in Sec. IV, effects of a phase transition is 



II. ASYMMETRIES IN COMOVER 
ABSORPTION 

We describe the evolution of the charmonium distribu- 
tion by the Boltzmann equation 
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where rii(ri,Pi,t) are the particle distribution functions, 
a^c is the comover absorption cross section for scatterings 
of charmonia tp with a comover c and v^ c their relative 
velocities. 

Let us first consider peripheral nuclear collisions where 
the overlap zone and the produced particle densities are 
small. Collisions are then rare and their effects can be 
calculated to first order with the free streaming distri- 
bution functions in the collision term of the Boltzmann 
equation Assuming Bjorken scaling along the beam or z- 
axis, v z = z/t. The free streaming distribution can at 
later times be approximated by |Io| , |lT[ 

n(x,p) ~ ^-S(y - S±(r v(r - r )) . (2) 

Here, r = (x, y) is the transverse radius and v = (v x ,v y ) 
the transverse velocity of the scatterers or comovers. As 
usual, i] — (1/2) log((t + z)/(i — z)) is the space-time ra- 
pidity and t = \ft 2 — z 2 is the invariant time, tq is the 
formation time for comovers. 

The transverse particle density distribution, S±(x, y), 
is for convenience approximated by simple gaussians 
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The rms radii of the overlap zones R x and R y at colli- 
sion are determined from nuclear thickness functions [|11| . 
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The x-direction is chosen along the impact parameter b 
and (x, z) constitutes the reaction plane with y perpen- 
dicular to the reaction plane. Replacing the gaussian 
transverse particle distribution by a constant density in- 
side an ellipse with axes R x _ y in Eq. ([|) does not change 
results by much. The only relevant scale is v x<y T / 'R x , y 
which contains the direction of pj_ and the transverse 
anisotropy. 

The integration over transverse coordinates and rapid- 
ity is straight forward (see Jll[] for details) . The scatterer 
momenta can also be summed over and in the follow- 
ing we will by (..) indicate averaging over scatterer mo- 
menta p c . Only an integral over proper time is left. With 
Ni = dNi/dp±dy we obtain for the comover absorption 
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where we have defined the dimensionless quantity 
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and the cross section times a transverse particle density, 
i.e. an effective "opacity" , as 
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In the participant or "wounded nucleon" model, which 
seems to describe the gross properties of relativistic 
nuclear collisions reasonably, both the rapidity density 
dN/dy and the transverse area R x R y scales approxi- 
mately with Ej>, except for very peripheral collisions (see, 
e.g., pjlfl). Therefore the opacity a of Eq. @ is approx- 
imately constant for central and semi-central collisions 
but decreases for very peripheral collisions. The almost 
constant opacity is different from the comover model of 
H which assumes densities to scale linearly with trans- 
verse energy. 

The dependence on the azimuthal angle 4> between the 
reaction plane and the charmonium transverse velocity 
(v x ,v y )^ = (cos 4>, sin <fi) , (see Fig. 1) is coupled to the 
source deformation 
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Defining the mean radius R 2 = (R 2 + R 2 )/2, we obtain 
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Here, we have expanded assuming S(v^/(v^ c )) <C I, 
which is valid for charmonia with p± ~ 1 GeV. 

The deformation can be approximated from the full 
transverse extent of the initial overlap of two nuclei with 
radius R A colliding with impact parameter b (see Fig. 1). 
They are simply: R x = R A -b/2 and R y = y 1 R\ - b 2 /^ 



and the corresponding deformation is 6 = b/2R A - Taking 
the rms radii of the nuclear overlap zone weighted with 
longitudinal thicknesses gives almost the same and may 
be related to the elliptic flow and HBT radii |pT|| . 

For peripheral collisions, where the comover absorp- 
tion is small, we obtain by integrating Eq. (Q) from the 
comover formation time tq to infinity 
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where we have introduced the function (C = 0.557..) 
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With typical relative velocities and formation times we 
find {v^ c }to <C R except for the very peripheral collisions. 
^From Eq. (||) thus £ 1 for central and semi-central 
collisions. 

The result of Eq. (|) for comover absorption gives the 
average suppression (see Fig. 2) as well as a characteris- 
tic oscillation with the azimuthal angle between reaction 
plane and the transverse momentum of the charmonia 
(see also Fig. 3). It is important to observe that the 
average and the amplitude of the absorption probes dif- 
ferent space-times of the collision and therefore different 
physics: The average absorption occurs at early times be- 
cause densities decrease inversely with proper time which 
leads to the logarithmic dependence of the initial time To 
as is also found in the comover model of ||. Contrar- 
ily, the amplitude of the oscillation is more sensitive to 
later times of order t ~ 2Rlv^, due to the extra factor 



(t — to) that appears in (|4|) when expanding for small 
deformations as in Eq. (||). The physics behind this late 
time effect on the amplitude, is that the particle has to 
travel a distance ~ R before the deformation of the source 
leads to more or less absorption. Likewise the amplitudes 
decrease with the transverse velocity as iA because the 
slow charmonia do not have the time travel the distance 
R before the comovers have expanded on a time scale 
R/v c . 

For more central collisions the free streaming assump- 
tion breaks down as interactions between comovers drive 
the system towards thermal equilibrium and hydrody- 
namic expansion. Also the absorption is a significant 
fraction and the linear expansion of Eq. (|^) breaks down. 
For simplicity, we shall assume that the free stream- 
ing result of Eq. (Q) is approximately valid for the co- 
movers in central collisions as well. The heavy charmonia 
do not thermalize but rather break up when interacting 
and therefore the distribution function for the remaining 
charmonia is free streaming to a good approximation. By 
taking proper time derivative of Eq. (^) and treat as 
a decreasing function with time, we obtain by integrating 
over time that 
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where the latter expansion is valid for central and 
semi-central collisions where (v^ c )tq/R and S are small. 
N^,(tq) is the number of J/iff's at time To where comover 
absorption sets in. 

In the wounded nucleon model also the transverse en- 
ergy Et scales with the transverse area, i.e., R x R y — 
R 2 oc Et- Consequently, the J /iff suppression decreases 
with centrality as a power law, J /iff oc E T a ^ 2 . It dif- 
fers from the exponential fall-off of the comover model of 
|^| that results from assuming that number densities and 
thus effective opacities scales linearly with Et- 

The transverse momentum, rapidity and centrality de- 
pendence of the amplitude in Eq. ( |12| ) differs from the 
average suppression. The amplitude is increase propor- 
tionally to dN/dy and p\ but decreases as S almost lin- 
early with centrality, Et- The relative amplitudes for 
J/iff, x, etc. are proportional to the absorption cross 
section and should therefore scale with these when cor- 
rected for fccddown whereas the DY background should 
not oscillate. 

We note that the prefactor to cos 4> is exactly the ellip- 
tic flow parameter (2^) if the absorption cross section is 
replaced by the momentum transport cross section, a tr 
|fll|| . The elliptic flow of pions and proton ~ 0.1 for 
pions and protons could be described with atr — 10 mb. 
Therefore, the magnitude of the charmonium absorption 
amplitude is approximately ~ 2v-ia^ c / atr —5% depend- 
ing on rapidity and transverse momentum. 



III. GLAUBER ABSORPTION 

Glauber absorption of charmonium on nucleons is be- 
lieved to be more important at SPS energies than co- 
mover absorption. The dissociation cross sections on 
nucleons and secondaries depend on various model as- 
suptions. Typically, for nuclear absorption one takes 
ct^at ~3-5 mb and assumes that comovers scatter with 
cross sections varying from a^ c ~ 1 — 5mb The 
relative amount of Glauber versus comover absorption 
is, however, unresolved at present since it is also possible 
to describe the charmonium by Glauber absorption -at 
least for p+A and S+U reactions- taking ctyjv —7 mb 
j|. Furthermore, it has been pointed out that feed-down 
effects from iff' and x _ decays into 3 /iff || or formation 
times Ji^JtJ can play a major role. 

As we shall concentrate on azimuthal asymmetries, we 
choose the standard model for Glauber absorption 

Gab (b) = Af J d 3 r d 3 r'p A (r)p B (r')6 2 (h - r ± + r' x ) 

xexp(-L A (r)-L B {r')) , (13) 



where pa.b are the nuclear densities, La(t) = 
cr^jv J z °° dz A pA( r ±, za), and analogously for Lb{t'). The 
normalization constant M is chosen such that Gab (b) = 
1 when La = Lb = 0. For simplicity we use the Glauber 
absorption of J /iff with uj/^jy = 3.6 mb as in D and no 
formation time. 

The straight line geometry implied by the Glauber 
model for absorption does not contain any azimuthal 
asymmetry, i.e., it is independent of <\>. In the NA50 
experiment the J/iff's are collected at pseudorapidies 
i] = 2.8 — 4 corresponding to angles of 2-7 degrees with re- 
spect to the beam axis. The J/iff's thus traverses a trans- 
verse distance r± ~ Rv± cosh(?y) ~ 0.2 — 0.9 fm inside the 
nuclei and their absorption might be affected by the den- 
sity distribution variation transversely. However, the net 
absorption from opposite direction depends on the cur- 
vature of and not the gradient of the transverse density 
distribution. Therefore, the resulting azimuthal depen- 
dence of Glauber absorption is ~ {r±/R x ) 2 — (r±/R y ) 2 ~ 
S(r±/R) 2 . It amounts to a few percent at midrapidity 
but is negligible at the forward rapidities. 

Momentum broadening affects the initial charmo- 
nium production rates and is another possible mecha- 
nism which could introduce azimuthally dependence in 
Glauber absorption. Momentum broadening has been 
found in pA and BA collision at relativistic energies 
and depends on the number of NN collisions and there- 
fore also the nuclear thickness function. The momentum 
broadening could be sensitive to the varying density pro- 
file depending on the range of the interaction, r^. As 
above the resulting azimuthal dependence of Glauber ab- 
sorption is ~ (TijR x ) 2 — (ri/Ry) 2 . Furthermore, it has 
to be multiplied with the relative momentum broaden- 
ing. With Ti ~ 1 fm we estimate that this effect is also 
small. Since DY pairs should be affected by momen- 
tum broadening of the cascading nucleons like Glauber 
apsorption but not by comover absorption, it is possi- 
ble to discriminate between the azimuthal dependence of 
Glauber and comover absorption by measuring the am- 
plitude for Drell-Yan (DY) production. 

These rough estimates clearly deserve more study in 
detailed microscopic models which, however, is outside 
the scope of this work. We emphasize that the time scales 
and distances involved for the initial Glauber scatterings 
are very short. Causality limits particles to probe trans- 
verse distance scales shorter than ~ Rz/l, where 7 is 
the Lorentz contraction factor. Such short distances are 
to be compared to the transverse radii, R x ,y, over which 
densities vary azimuthally. 

Predictions for RHIC and LHC energies depend on 
the charmonium absorption cross section at these ener- 
gies. If formation times are t^^R/j, the Glauber ab- 
sorption might dissappear or be enhanced at higher en- 
ergies due to time dilation and different production 
mechanisms as e.g. suggested by the color-singlet or the 
color-octet model. However, pA Fermilab data on J /iff 
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production |13|] indicates that the absorption cross sec- 
tion does not change going from 200 to 800 GeV lab 
energy. For the estimates presented below we therefore 
assume that the Glauber absorption is unchanged going 
from SPS to RHIC and LHC energies. Furthermore, the 
minor azimuthal dependence of Glauber absorption will 
be ignored. 



IV. ESTIMATES FOR SPS, RHIC AND LHC 

First we estimate the parameters in play. The un- 
known ones are the comover absorption cross section, 
tr^c, and the formation time To. The other parameters 
and their dependence on impact parameter can be es- 
timated from experiments and results in an almost con- 
stant opacity as explained above. We give an estimate for 
central Pb + Pb collisions at SPS energies. Here, the rms 
radius is R 2 ~ 10fm 2 , and the total rapidity density is 
dN tot jdy ~ 500. One can assume [|J that only heavy res- 
onances as p,uj,(j), etc. have sufficient energy to destroy 
the charmonia and take dN c /dy ~ 0.5dN tot /dy. The 
average relative velocity for charmonia with transverse 
momentum ~ 1 GeV/c is v^ c ~ 0.6c. A comover absorp- 
tion cross section of order tr^ c = (2/i)a^^ ~ 2.4 mb is 
typically assumed, and so we obtain the opacity a ~ 0.3 
for central and semi-central collisions. Its dependence 
on centrality for peripheral collisions is calculated in the 
wounded nucleon model 

At RHIC and LHC energies the rapidity density is ex- 
pected to increase with collision energy and we shall sim- 
ply assume that they double at RHIC energies and triple 
at LHC energies. Consequently, the opacity is a =0.3, 
0.6, and 0.9 at SPS, RHIC and LHC energies respectively. 
In Fig. 2. we show the cor resp onding comover suppres- 
sion of the J/4" from Eq. (|TlP^2[) assuming To = lfm/c 
and taking the same standard Glauber absorption from 
Eq. ( |l3] ) with (j^at = 3.6 mb at all three energies. 

The modulation of the .//^-suppression with angle be- 
tween pj_ and reaction plane is shown in Fig. 3 for 
five centralities or Et corresponding to impact param- 
eter b/R ~ 2, 1.5, 1, 0.5, 0. For a prediction at RHIC we 
take the same Glauber absorption as in Fig. (2) as well 
as the SPS values for v$ ~ 0.3c and v^, c ~ 0.6c. How- 
ever, a ~ 0.6 is employed for RHIC energies due to the 
larger rapidity densities. The relative amplitude is given 
by Eq. (|l^) and decreases with centrality as the defor- 
mation. For near-central collisions the source expands 
significantly and the deformation may be smaller than 
that at initial overlap, 6 ~ b/2R. The magnitude of the 
amplitude is up to ±7% for semi-peripheral events. As 
the relative amplitude scales with the rapidity density, it 
is smaller by a factor of two at SPS energies but larger 
at LHC energies. 

The relative amplitudes should be measurable by tech- 
niques similar to those used for extracting elliptic flow M . 
By taking the ratio of the amplitude to the average JpS> 



or using a minimum bias analysis, the largest statistical 
uncertainty (the DY background) is removed. The NA50 
minimum bias data has a relative error bar of 1-2%, i.e., 
an order of magnitude less than the variation in the J "/ W 
absorption, [N$(<f> = 0) — NjAcf) = 7r/2)]/A^, which is 
twice the amplitude in Eq. (|12[). 

It would also be most interesting to measure and com- 
pare the amplitude for the various charmonium states 
J/ty, x as we U as f° r the DY background. The 'I'' is 
expected to have a larger dissociation cross section and 
the amplitude should be correspondingly larger. The am- 
plitude for the DY background should vanish as discussed 
in the previous section. 



V. EFFECTS OF A PHASE TRANSITION 

The anomalous J/^ suppression measured by NA50 is 
possible due to the formation of a plasma state leading to 
color screening of cc pairs [^]|| . As the plasma is created 
at high energy densities present at early times it would 
increase the average suppression. However, the ampli- 
tude would be less affected since it is caused by late time 
absorption. We can mimic this situation by employing a 
larger absorption cross section at early times than at late 
times or in terms of Eq. (^) a larger a in the average 
power law suppression than in the prefactor to the am- 
plitude. Both quantities can be measured as function of 
centrality and it would be most interesting to follow their 
behavior for the J/^f and suppression as well as the 
subthreshold enhancement - in particular in the regions 
of anomalous J /iff and W and suppression. 

If the charmonium absorption cross sections are 
strongly energy dependent 0], the average cross section 
(ifycCi/ic) may also be larger at earlier times since the rela- 
tive energies are larger here. In the above calculations we 
ignored relative longitudinal velocities by assuming the 
8{y — rj) factor in the distribution. This approximation 
breaks down at early times and one should therefore take 
an effective cross section at larger collision energy in Eq. 
(|l2| ) at early times. However, whereas the energy depen- 
dent cross sections should lead to a smooth variation of 
the average J/^-suppression as well as its amplitude, the 
plasma leads to an anomalous and sudden suppression at 
a certain centrality or Et f|||]- 

VI. SUMMARY 

In summary, we have investigated suppression of char- 
monia in non-central collisions and find a characteristic 
oscillation with azimuthal angle between transverse mo- 
mentum and reaction plane due to comover absorption. 
Glauber absorption was estimated to contribute less to 
the amplitude because causality strongly constrains the 
transverse distances probed and the resulting azimuthal 
dependence. Therefore, the azimuthal variation may be 
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exploited to gain more information on the suppression 
mechanisms. 

Estimates for comover absorption predicts an az- 
imuthal variation of order 5-20% in the number of J/^f 
depending on centrality, size and energy of the colliding 
nuclei. Such an oscillation should be detectable (if the 
reaction plane can be determined event-by-event) as it 
is not sensitive to the statistical uncertainty of the DY 
background. In comparison, the current minimum bias 
data of NA50 has error bars of 1-2%. 

The amplitude has a strong dependence on rapidity 
and transverse momentum of the charmonia and on cen- 
trality that is different from the average suppression. 
The amplitude is proportional to dN/dy and p\ and de- 
creases almost linearly with centrality Et- The relative 
amplitudes for J/ty, ip' and \ are proportional to the 
absorption cross section and should therefore scale with 
these when corrected for feeddown. The DY background 
should, however, not oscillate. 

The average absorption and its amplitude are more 
sensitive to early and later times respectively. This dis- 
tinction may further be exploited to obtain additional 
information on the absorption mechanism and whether 
an anomalous suppression is due to a phase transition. 
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FIG. 1. Reaction plane of a semi-central Pb + Pb collision with impact parameter 6 = Rpb — 8fm. The overlap zone is 
deformed with R x < R y . The reaction plane (x, z) is rotated by the angle (p with respect to the (outward) transverse particle 
momentum px- The scatter plot shows charmonium dissociation points in a microscopic cascade simulation Jl3] at SPS energies. 




FIG. 2. Comover and Glauber suppression of J/*, Eqs. (0-0), shown by curves as function of transverse energy Et- The 
opacities are: a = 0.0, 0.3, 0.6, 0.9 corresponding approximately to pure Glauber absorption, SPS, RHIC and LHC energies (see 
text for details). Also shown is NA50 data |l[ normalized to DY and minimum bias (open and filled circles respectively). 
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FIG. 3. J /'if suppression as function of azimuthal angle cf> normalized to its average for various centralities from Eq. (|l2"|). 
The RHIC parameters a = 0.6 are employed (see text). 
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